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Abstract 

We report on measurements of the inclusive production rate of S + and S° 
baryons in hadronic Z decays collected with the L3 detector at LEP. The S + baryons 
are detected through the decay S + — > p7r°, while the S° baryons are detected via 
the decay mode S° — > A7. The average numbers of S + and S° per hadronic Z decay 
are measured to be: 

(N s+ ) + (N t+ ) = 0.114 ±0.011 stai ± 0.009 5 ^ 
(JVso) + (iVgo) = 0.095 ±0M5 stat ±0M3 syst . 

These rates are found to be higher than the predictions from Monte Carlo hadroniza- 
tion models and analytical parameterizations of strange baryon production. 
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Introduction 



Measurements of hadron production in e + e~ annihilation are important to understand the frag- 
mentation process of quarks and gluons into hadrons. Perturbation theory cannot be applied 
to describe this process and the theoretical description is based on phenomenological models 
as implemented, for example, in the JETSET 0, HERWIG [2 or ARIADNE Monte Carlo 



Recently, two analytical models have been proposed to describe the hadron production rates 
in e + e~ annihilation [HE]. The model of Reference [3], referred to from here on as the "string- 
based model," is derived partially from string fragmentation and describes the production of 
light mesons and baryons using the simple formula: 



where (N) is the number of primary hadrons of spin J, containing N s strange quarks, produced 
directly by string fragmentation. The quantity -Ebind is the hadron binding energy: -Ebind = 
Miadron — J2i m qi > where m qi are constituent quark masses, and T is the effective temperature 
of hadronization in this model. The overall normalization of the rate is controlled by C. 
The evolution of C with center-of-mass energy is predicted by perturbative QCD j^J. The 
parameter j s is a universal suppression factor for strange quark production, while is a 
suppression factor of baryon production relative to meson production. This model has five free 
parameters: C at some fixed energy scale, Cb, Is, T and Am = m s — m u , the mass difference 
between strange and up quarks. The model of Reference [Sj, referred to from here on as the 
"hadron gas model," is based on thermodynamical considerations and has three adjustable 
parameters: the total volume V of the hadron gas, a strangeness suppression factor similar 
to 7<- in Reference j3], and the hadronization temperature T. In the hadron gas model, no 
additional parameters are needed to reproduce the relative production rate of baryons versus 
mesons. Both of these models give a good description of the full spectrum of light hadron 
production in addition to predicting the rates for strange baryons. The string-based model also 
describes well the relative production rate of hadrons containing b and c quarks. In order to 
compare the experimental data to these two analytical models which predict primary hadron 
production, the observed rate must first be corrected for two effects: (i) hadrons containing 
quarks produced in the primary electroweak interaction, (ii) 'feed-down' due to the decays of 
unstable hadrons produced in string fragmentation [4J or, correspondingly, the phase change 
that produces the hadron gas [3]. The feed-down correction is made using known hadron 
branching ratios [7j. 

This paper presents measurements of S + and S° strange baryons produced in hadronic 
Z decays at center-of-mass energies in the range y/s = 91 ± 2 GeV, from 71.1 pb _1 of data 
collected by the L3 detector during the LEP running periods of 1994 and 1995. The S + and S° 
baryons provide a good test of the hadron rate models because of their relatively low feed-down 
correction as compared to the lighter p and A baryons. The use of the decay modes which have 
photons, S + — > p7T° and S° — > A7, enable the measurements to be made down to low baryon 
momenta due to the low-energy photon detection capability of the L3 detector. 



generators. 
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The L3 Detector 

A detailed description of the L3 detector and its performance is given in Reference [S]. This 
analysis makes use of the electromagnetic calorimeter made of Bismuth Germanate (BGO) crys- 
tals, which accurately identifies and measures photons from an energy of 50 MeV to 200 GeV. 
In order to benefit from the low-energy photon detection of the calorimeter, its performance 
was studied and parameterized The resulting resolution functions are shown in Figure ^ 
In addition, corrections due to electronic noise and calibration accuracy are taken into account. 
The electronic noise terms are yN • cr intr with cr intr = 1 MeV and iV ■ a corr with a corr = 1.6 MeV 
where a[ ntT is the mean intrinsic noise per channel and <r corr is the mean correlated noise per 
channel, and N is the number of crystals used to compute the shower energy. The calibration 
accuracy term for 1994 and 1995 is a C3 \ ■ E with <r ca i = 1.5% as estimated from Bhabha events. 
These terms are added in quadrature to the resolution functions of Figure [T] 

Hadronic Event Selection 

The inclusive production rate measurements are based on a selection of hadronic Z decays 
which have a measured primary vertex. A precise measurement of the primary vertex in all 
three dimensions is achieved in more than 99% of all hadronic events. The hadronic selection 
requires a well-balanced, high particle-multiplicity event with the full collision energy measured 
by the detector and relies on the following criteria: 

1. The total energy observed in the detector, including the momenta of muons measured in 
the muon spectrometer, E w - lS) is restricted to the range 0.6 < E vis / y/s < 1.4. 

2. The energy imbalance along the beam direction, E\\, must satisfy \E\\ \ / E vis < 0.4. 

3. The transverse energy imbalance, E±, with respect to the beam direction, must satisfy 
E ± /E vis < 0.4. 

4. The number of energy clusters reconstructed in the calorimeters is required to be larger 
than 12. 

A data set of 1.6 million hadronic Z decays was selected for this analysis. 

Analysis Procedure 

The events passing the hadronic event selection are analyzed for the inclusive production of 
E + and S° in the decay modes p7r°(7r° — > 77) and A(— ► p7r~)7, respectively. The p7r° and A7 
mass distributions are computed, and a fitting procedure is applied to count the number of 
selected candidates in the mass peak. The anti-particles of the S + and S° are included in the 
corresponding distributions to increase the statistical significance of the data samples. 

Central to the selections of the S + and S° candidates are the kinematical constraints de- 
scribing their decay. A probability is assigned to each candidate based on the minimum \ 2 
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which satisfies the decay constraints. The x 2 has the form: 

-^tracks 

X 2 = Yl (** ~ f^s , <?i ) ) 1 (ti - f(v 8 , & ) ) + (V - V cvcn t ) T Kvent ~ *Wt ) 

t=l 

-^clusters 

+ («i - <f(i£, q^)) T E~ l {ui - g(v s , g)) (2) 

i=i 

where the number of charged-particle tracks and electromagnetic clusters used to reconstruct 
the decay are respectively iV trac ks and iV c i us ters; U is the vector of measured parameters for the 
2 th track and Gi the corresponding covariance matrix; f{v S) qi) are the fit parameters assuming 
the track originated from the secondary vertex v s with momentum qf, tWnt is the event- vertex 
with covariance matrix V; v is the fit to the production vertex; Ui is the vector of measured 
parameters for the i th electromagnetic cluster with covariance matrix Ef, and g{vl,qi) are the 
fit parameters assuming that a photon originated from the secondary vertex v* s with momentum 
qi. The details of the constrained fitting technique are found in Reference In addition to 
the constraint probability cuts, other selection criteria are also applied to reduce combinatorial 
background. 

The fitting procedure of the invariant mass distributions consists of several steps. The 
JETSET Monte Carlo simulation is used to predict the shape of the p7r° and A7 mass distri- 
butions separately for the background and the S + and S° signals. The Monte Carlo signal and 
background distributions are binned into histograms which are then smoothed with a spline 
fit. The spline fit reduces background fluctuations in the signal region by interpolating between 
the sidebands of the signal. The use of binned histograms with bin widths comparable to the 
expected mass resolution of the E + and E° signals is found to be optimal for exploiting the 
histogram shape differences between signal and background. A maximum likelihood fit is per- 
formed for both the signal and the background normalizations, simultaneously. The likelihood 
is constructed by computing the Poisson probability for each bin and then making the product 
over all bins in the histogram, as follows, 

^bin s / f , \di 

L{s , b) = e ^n^0_ (3) 

where s and b are respectively the total number of signal and background events, is the 
number of observed data events in a given bin and fi and gi are the number of expected events 
in a given bin for signal and background, respectively. 

The production rates of the S + and S° are determined by a kinematic extrapolation which 
assumes the JETSET Monte Carlo production distributions. 



H + Selection 

The signature of a S + — > p7r°(7r° — > 77) decay is a single charged track and two photons 
measured in the calorimeter as shown in Figure The kinematical information of the decay 
configuration, in addition to the assumption that the S + is produced at the primary vertex, 
can be expressed as follows: 

d' x [p p + p 7l + p 72 ] = (4) 
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m 77 = m n o (5) 

where the vector d' is the direction of flight of the S + at the time of decay; the vectors p p , p 7l 
and p 72 are respectively the proton momenta and the momenta of the two photons; and m 77 and 
m n o are the calculated mass of the diphoton system and the 7r° mass, respectively. A probability, 
V sv tx, is assigned to each S + candidate based on the minimum \ 2 which satisfies constraints (J1J) 
and (jHJ). The x 2 is obtained from equation (J2J). The computation of the constraint probability 
requires the accurate description of the photon resolution in energy and angle in the calorimeter. 
The selection of S + candidates consists of the following cuts: 

1. The constraint probability, V sv t x , must be larger than 0.2. 

2. The energies of the two photons must be greater than 55 MeV and their polar angle must 
be in the BGO barrel acceptance (|cos# 7 | < 0.74). 

3. The 7T° boost, 7^0 = E^o/m^o, is required to be larger than 2.0. 

4. The transverse momentum of the proton must exceed 800 MeV. 

5. The transverse decay radius, r_i_, of the S + is required to be in the range 5 mm < r± < 
70 mm ensuring the decay to be inside the inner radius of the Silicon Microvertex Detector 
(SMD). 

6. Charged tracks which are identified as coming from A or K° decay are rejected as proton 
candidates from S + decay. 

Figure El shows the mass distribution of the selected p7r° candidates. The number of selected 
S + — > p7r° decays determined by the fit is found to be 342 ± 33. 



E Selection 

The S° is measured in the decay mode S° — > A7. This mode accounts for nearly 100% of 
the E° branching fractions. Identification of the A is done using the p7r~ decay mode. The 
proton is assumed to be the track with the highest momentum of the two tracks forming the A 
candidate. The kinematic constraints of the S° — > A(— >• p7r~)7 decay, shown in Figure Eb, are 
summarized by the following equations: 

d x [p p + pV] = (6) 

m pn - =m A (7) 

where the vector d points from the primary vertex to the A decay point; the vectors p p and 
p„— are respectively the proton and pion momenta; and m pn - and are the calculated mass 
of the proton- pion system and the A mass, respectively. A probability, Pi v tx, is assigned to each 
A candidate based on the minimum \ 2 (equation EJ) which satisfies constraints (0) and (J7j). A 
large combinatorial background to the A signal is present in the limit that 

v s = v event (8) 

where v s is the A decay vertex and invent is the primary event vertex. A constraint probabil- 
ity, "Ppvtx, is also computed for A candidates to satisfy equation (jHJ). The A signal has less 
background when V pvtx is small. 

The A selection consists of the following cuts: 
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1. The transverse momenta of the proton and pion coming from the A decay must both 
exceed 200 MeV. 

2. The constraint probability, Pi v t x > must be larger than 0.01, and the constraint probability, 
P pvtx , must be smaller than 10~ 12 . 

3. The transverse decay radius, r±, of the A is limited to the range r± < 70 mm ensuring 
the decay to be inside the inner radius of the SMD. 

After applying the A selection cuts with the exception of constraint (JZj), the reconstructed 
p7r~ mass spectrum is plotted in Figure |U where the A is clearly visible. To reconstruct the S° 
decay, an electromagnetic cluster in the BGO detector is combined with a selected A. 

The following additional criteria is applied to select S° candidates: 

1. The energy of the photon, E 11 is restricted to the range 55 MeV < E 1 < 1 GeV and the 
photon polar angle must be within the BGO barrel acceptance (|cos# 7 | < 0.74). 

2. The angle, 9\, of the A with respect to the S° momentum in the S° rest frame is required 
to satisfy cos 9\ < 0.2. 

Figure El shows the mass difference distribution of the selected A7 candidates. The number 
of selected S° — > A(— > pvr~)7 decays is found to be 263 ± 42. 

Rate Measurements 

A total of 6.2 million hadronic Z decays were simulated with the JETSET Monte Carlo to 
estimate detector effects. The simulation accounts for the run conditions in luminosity-weighted 
periods throughout the data-taking. 

For signal and background fitting, further tuning of Monte Carlo distributions is performed. 
The signal shape predicted by the simulation is adjusted to fit the observed peak position and 
resolution measured in the data. The change in the rate for a corresponding change of 0.5 in 
the negative log likelihood in the fit of the background and signal normalizations is taken as a 
systematic error. 

To estimate the statistical error on the Monte Carlo determination of the background shape, 
the signal is fit with a fixed normalization of the background corresponding to two cases: 
a) increasing the background-level by one standard deviation (given by the two parameter 
fit), and b) decreasing the background-level by the same amount. In the case of the S°, 
where the photon energy spectrum goes down to threshold, the observed difference in the low- 
energy photon energy spectrum between Monte Carlo and data is used to reweight the mass 
difference distribution. The correlation on the rate measurement between reweighting based on 
background versus background plus signal is quoted as a systematic error. 

In order to estimate the error on the kinematic extrapolation of the observed rate to the full 
phase space of the S + and E° baryons produced in the fragmentation process, several Monte 
Carlo generators are used to predict the extrapolation The variations in the fractions 

of decays inside the kinematic ranges of the E + and S° selections are given in Table d These 
are taken as systematic errors. Varying the S + and S° selection cuts over a wide range of cut 
values has no observed systematic effect on the measured rates. 

Using the JETSET Monte Carlo as a reference for the extrapolation, the efficiencies for 
S + and S° are found to be 0.19% and 0.18%, respectively. The error on the efficiencies from 
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Monte Carlo statistics is taken as a systematic error. The sources of systematic errors in the 
production rate measurements are listed in Table El The average numbers of S + and S° per 
hadronic Z decay are measured to be: 



Our measurements of the S + and S° production rates in hadronic Z decays are compared 
to baryon production rates predicted by JETSET, HERWIG and ARIADNE models in Ta- 
ble El PUJ. In this table, the numbers obtained from the string-based and the hadron gas 
models are also listed. The predictions of all of these models underestimate our measured 
rates. Our measurements are consistent with other measurements performed at LEP [TlTll3| . 
However, we observe rates which are somewhat higher. This difference could arise from the 
low momentum part of the baryon production spectra. The use of low energy photons in the 
measurements of the S + and S° production rates allows the kinematical ranges covered by the 
L3 measurements to extend down to lower baryon momenta than those measured by the other 
LEP detectors. 
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A/Tonfp Hariri 


f inside 


f inside 

is 


ARTADNF, 4 08 


o 31 7 


100 


JETSET 7 4 


0.321 


0.097 


HERWIG 5.9 


0.311 


0.091 


Maximum Spread 


3.1% 


9.0% 



Table 1: Monte Carlo predictions of the fractions of decays inside the kinematic ranges of the 
E + and S° selections. 



Source of 
Systematic Error 


(%) 


a (Zh(^o) (%) 


Signal Shape 
Resolution 


6.1 


7.0 


signal oliape 
Peak Position 


0.6 


1.9 


Background Shape 
Statistics 


2.4 


6.4 


Background 
Reweighting Procedure 




2.9 


Kinematic 
Extrapolation 


3.1 


9.0 


Monte Carlo 
Signal Efficiency 


3.5 


3.8 


Total 
Systematic Error 


8.1 


14.1 



Table 2: Sources of systematic error in the production rate measurements. 
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Monte Carlo 


(7V S+ ) + (7V S+ ) 


(7V s o) + (iVgo) 


JETSET 


0.0711 


0.0691 


HERWIG 


0.0978 


0.0692 


ARIADNE 


0.0716 


0.0685 


string-based model 


0.0740 


0.0774 


hadron gas model 


0.0732 


0.0766 



Table 3: Model predictions of S + and S° production rates in e + e~ annihilations at 91 GeV. 
The predictions of the JETSET, HERWIG and ARIADNE models have been obtained with 
parameters tuned to L3 global event shape distributions and the charged-particle multiplicity 
distribution. 
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Figure 1: Resolution of electromagnetic showers in energy and azimuthal angle <fi. In plot a) 
the two curves describe the energy resolution in the BGO barrel. The full resolution curve 
accounts for the fraction of measured photons which lack full shower containment. These give 
a tailing distribution with less resolution than the intrinsic Gaussian resolution curve. In plot 
b) is the 0-resolution function. At low-energy, the angular resolution is limited by the crystal 
size. The dotted line shows the extrapolated angular resolution when this effect is ignored. 
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Vertex <=s 

Figure 2: The decay of a S + is shown in plot a). The position of the decay vertex, v S) and the 
path of the S + from the primary vertex to the decay vertex are indicated. The vector d' is the 
direction of flight of the S + at the time of decay. The decay of a S° is shown in plot b). The 
point v s indicates the position of the decay vertex of the A, and the vector d is the direction of 
flight of the A. 
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Figure 3: The mass distribution of S + candidates. The shapes of the signal and of the back- 
ground are obtained from Monte Carlo simulation as described in the text. 
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Figure 4: The p7r~ mass distribution showing the A candidates used to measure the S°. 
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Figure 5: The mass difference distribution M(A7) — M(A) of S° candidates. The shapes of the 
signal and of the background are obtained from Monte Carlo simulation as described in the 
text. 
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